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We investigate the existence and the properties of fully separable (fully factorized) ground states 
in quantum spin systems. Exploiting techniques of quantum information and entanglement theory 
we extend a recently introduced method and construct a general, self-contained theory of ground 
state factorization in frustration free quantum spin models defined on lattices in any spatial dimen- 
sion and for interactions of arbitrary range. We show that, quite generally, non exactly solvable 
translationally invariant models in presence of an external uniform magnetic field can admit exact, 
fully factorized ground state solutions. Unentangled ground states occur at finite values of the 
Hamiltonian parameters satisfying well defined balancing conditions between the applied field and 
the interaction strengths. These conditions are analytically determined together with the type of 
magnetic orderings compatible with factorization and the corresponding values of the fundamental 
observables such as energy and magnetization. The method is applied to a series of examples of in- 
creasing complexity, including translationally-invariant models with short, long, and infinite ranges 
of interaction, as well as systems with spatial anisotropics, in low and higher dimensions. We also 
illustrate how the general method, besides yielding a large series of novel exact results for complex 
models in any dimension, recovers, as particular cases, the results previously achieved on simple 
models in low dimensions exploiting direct methods based on factorized mean-field ansatz. 

PACS numbers: 75.10.Jm, 03.65.Ca, 03.67.Mn, 64.70.Tg 



I. INTRODUCTION 

Quantum information theory is an area of scientific 
investigation that has witnessed an enormous progress 
in the last decade [ij. In the framework of quantum 
information science, paradigmatic systems of condensed 
matter physics such as, e.g., spin chains and harmonic 
lattices, are analyzed from the point of view of their in- 
formation content and ability to process and transfer in- 
formation. Fundamental concepts of statistical mechan- 
ics and probability, such as the Shannon-von Neumann 
entropy 0], play a central role in the quantification of 
bipartite quantum entanglement 0]. On the other hand, 
exciting advances in quantum information research are 
in a sense paying back the debt: The mathematical and 
theoretical toolkit developed for the characterization and 
quantification of quantum entanglement has proven use- 
ful to tackle questions and improve our understanding 
in the investigation of strongly correlated systems and 
quantum phase transitions [J]. 

Perhaps the most interesting development stemming 
from the application of the tools of entanglement the- 
ory concerns the study of correlation scaling in diversely 
connected systems, with its relations to conformal field 
theory and the proof of longstanding conjectures on en- 
tropic area laws in spin and harmonic lattice systems Q . 
These investigations and the associated results are of fun- 
damental interest because they allow to enhance the con- 
trol and establish firmer bounds and limits of applicabil- 



ity on many important tools of simulation and numerical 
analysis of complex many-body systems, a fundamental 
task in condensed matter, such as the density matrix 
renormalization group (DMRG) algorithms [|, matrix 
product representations (3], the entanglement renormal- 
ization ansatz and weighted graph states methods 
i- 

The pioneering application of quantum information 
concepts to condensed matter was the observation that 
two-body entanglement, as quantified by the concur- 
rence, in the ground state of a cooperative system, ex- 
hibits peculiar scaling features approaching a quantum 
critical point [l^. These seminal studies helped to clar- 
ify that at quantum phase transitions [llj, the dramatic 
change in the ground state of a many-body system is 
associated to, or reflected by, a change in the way that 
quantum fluctuations are correlated, that is, in the way 
entanglement is distributed among the elementary con- 
stituents. Many recent efforts have thus been aimed at 
understanding the behavior of different measures of en- 
tanglement at quantum phase transitions, and assessing 
the corresponding enhancement of properties useful for 
applications in quantum technology and quantum engi- 
neering m. However, since most entanglement measures 
can be often rewritten in terms of the conventional n- 
point correlation functions, the presence of a divergence 
in the former at quantum criticality can be directly traced 
back to a divergence in the latter although in some 
particular instances entanglement-based studies allowed 
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to discover novel types of phase transitions [TjI ■ 

Finally, a more technologically oriented product of the 
interaction between these two areas of research is the 
study of quantum spin chains as natural information car- 
riers and distributors [14 L able to realize tasks such as 
quantum state transfer and storage possibly 
with maximal fidelity in chains with suitably engineered 
couplings and dynamics. 

Quantum engineering and quantum technology are be- 
ing indeed developed at a fast pace. Quantum devices are 
vigorously pursued for applications ranging from nano- 
sciences to quantum computation and entanglement- 
enhanced metrology . Despite a large variety of possi- 
ble implementations involving different physical systems, 
many relevant properties of such devices can be in fact 
investigated in a unified setting by appropriate mappings 
to quantum spin models [l^, [l^ . Control of separabil- 
ity and entanglement in ground states of quantum spin 
models plays an important role in entanglement-enabled 
quantum technology applications, and needs to be char- 
acterized and quantified in different physical regimes [J. 
From a more fundamental perspective, the determination 
of exact solutions endowed with precisely known prop- 
erties of separability or entanglement, can be of great 
relevance in the study of advanced models of condensed 
matter and cooperative systems that are in general not 
exactly solvable. 

In this paper we exploit quantum information tools to 
tackle another important problem in condensed matter 
physics: the determination and characterization of exact 
ground states of cooperative systems characterized by the 
property of being in the form of a fully factorized, tensor 
product of single-particle pure states, with no quantum 
correlations between the individual constituents. Histor- 
ically, the occurrence of totally factorized (unentangled) 
ground states of quantum many-body systems was first 
discovered in the one-dimensional anisotropic Heisenberg 
model with nearest-neighbor interactions by Kurmann, 
Thomas and Miiller [20|, by adopting a direct method 
in terms of product state ansatz. This result was later 
re-derived and extended to two dimensions using quan- 
tum Monte Carlo numerical methods (2l|. The direct 



product-state ansatz method of Kurmann, Thomas and 
Miiller has been later extended to anisotropic rings with 
interactions of arbitrary range, by exploiting the observa- 
tion that factorized ground states break the parity sym- 
metry j23j ■ Complex quantum systems exhibiting 
cooperative behaviors, whose ground states are typically 
entangled [2^ . may thus admit, for some non trivial val- 
ues of the Hamiltonian parameters, a ground state which 
is completely separable. The occurrence of factorization 
at finite or even strong values of the couplings is thus an 
effect of a delicate balancing between interactions and 
external fields. 

The phenomenon of ground state factorization is par- 
ticularly intriguing as it appears to be associated with the 
presence of an "entanglement phase transition" with no 
classical counterpart [25|. Furthermore, for the purposes 



of quantum engineering applications, that employ dis- 
tributed entanglement in order to manipulate and trans- 
fer information p^ . factorization points need to be ex- 
actly identified and either avoided, in order to guaran- 
tee the reliable implementation of quantum devices, or 
properly exploited for the dynamical creation of strongly 
entangled multipartite states of large assemblies of mi- 
crosystems (graph and cluster states). Finally, and also 
quite importantly, for models not admitting exact general 
solutions, achieving knowledge of the exact ground state, 
even if only for the restricted nontrivial set of parameters 
associated to factorization, would allow (i) to prove the 
existence of an ordered phase and characterize it; (ii) to 
build variational or perturbative approximation schemes 
around the exact factorized solution, that may then be 
used as test benchmarks for the validity and the precision 
of numerical algorithms and simulations. Unfortunately, 
the direct method based on the ansatz of a factorized 
solution is neither sufRciently general nor mathemati- 
cally tractable, apart from the notable exceptions of one- 
dimensional spin chains with short-range (20| . infinite- 
range |26|. and some types of long-range exchange inter- 
actions [ll, [H . 

In a recent work [27t we introduced an analytic method 
that allows to determine exactly the existence of factor- 
ized ground states and to characterize their properties in 
quantum spin models defined on regular lattices of any 
size (finite as well as at the thermodynamic limit) , in any 
spatial dimension, and with spin-spin interactions of arbi- 
trary range. In correspondence to rigorously established 
ground state factorizability, the method also allows to 
determine novel sets of exact solutions in generally non 
exactly solvable models. All the previous particular find- 
ings can be rigorously re-derived and extended within a 
unified framework inspired by concepts of quantum in- 
formation science. 

In the present work we extend the original method 
introduced in Ref. and we develop a general 

self-contained theory of ground state factorization in 
(frustration-free) quantum spin systems: With no a pri- 
ori assumption on the nature of the magnetic ordering, 
we derive a complete and closed set of conditions that 
have to be satisfied by the Hamiltonian interaction pa- 
rameters in order for the ground state to be factorized, 
at precisely determined values of the external magnetic 
field, with well defined values of the magnetization and of 
the energy. Besides the completely rigorous re-derivation 
of the few previously known results on quantum factor- 
ization points [13, [21I, [13, [H, [13], our general analytic 
method allows to determine novel, exact factorization 
points and factorized ground states in, generally non ex- 
actly solvable, spin-1/2 models as well as in the corre- 
sponding Hamiltonian models of any higher spin. These 
results hold true regardless of spatial dimensionality and 
interaction range, and can be easily extended even to 
non-translationally invariant systems. We illustrate the 
versatility of the method with a series of different appli- 
cations, including models with nearest-neighbor interac- 
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tions (for which wc provide the full factorization diagram, 
generalizing the results of Refs. [13, [H, HI] and of our 
recent work [l^l), models with long-range interactions 
(notably including, among others, the fully connected 
infinite-range Lipkin-Meshkov-Glick (LMG) model [1^ 
and systems with unbalanced next-to-nearest-neighbor 
interactions on cubic lattices), and models with spatial 
anisotropics. 

The paper is organized as follows. An excursus on 
the problem of ground state factorization is presented in 
Sec. ini In Sec. IIIII we recall the quantum information 
toolkit upon which our method is based, in particular 
the formalism of single-spin, or single-qubit, unitary op- 
erations and associated entanglement excitation energies, 
previously introduced for the characterization and quan- 
tification of bipartite entanglement from an abstract ge- 
ometrical perspective [1^ [23 • The general method for 
the rigorous determination of factorized ground states, 
that generalizes the scheme introduced in Ref . [l^] , is dis- 
cussed in full detail in Sec.lIVI Relevant applications and 
examples, including models defined on lattices of different 
spatial dimensions and with different interaction ranges 
as well as extensions to models with spatial anisotropics 
and with dimerized interactions are discussed in Sec. fVl 
SeclVH Scc. lVIIl and Sec. lVIIll Concluding remarks and 
outlooks on future perspectives are presented in Sec. IIXI 



II. THEORY OF GROUND STATE 
FACTORIZATION: STATE OF THE ART 

Quantum fluctuations in the ground states of cooper- 
ative many-body systems are typically highly correlated, 
that is, those states contain in the most general instance 
a strong degree of entanglement distributed among the 
individual components. This is one of the reasons why 
their exact characterization is often hard to be accom- 
plished by analytical or even numerical methods. How- 
ever, there may exist particular solutions corresponding 
to special values of the Hamiltonian parameters such 
that the ground state becomes "classical-like", in the 
sense of being a full product state of single-site factors 
(completely separable, fully factorized state). It is quite 
surprising that despite the inherent appeal of such ex- 
act, special solutions, and the renewed interest stemming 
from possible applications to quantum technology, very 
few analytical advances were obtained in the last quarter 
of a century. 



A. The direct method: Product state ansatz 

The first systematic work on ground state factoriza- 
tion in quantum spin models was completed by Kur- 
mann, Thomas, and Miiller in 1982 [2^. They focused 
on the one-dimensional nearest-neighbor anisotropic 
Heisenberg-like antiferromagnetic spin- 1/2 model in 
transverse field, with the assumption that all the spin- 



spin couplings take nonnegative values. Kurmann et al. 
identified a precise value of the external magnetic field, 
the "factorizing field" , corresponding to which the system 
admits a fully factorized ground state. They employed 
a very direct, "brute force" method, i.e. formulating a 
product state ansatz and verifying if and for what values 
of the Hamiltonian parameters it satisfies the stationary 
Schrodingcr problem with lowest energy. Even though 
this approach appears obvious and extremely easy to pur- 
sue, the positive result obtained by Kurmann et al. in the 
case of the 1-D Heisenberg model was enabled by some 
special properties of the case they analyzed. In fact, in 
order to verify the eigenvalue equation for the product 
state ansatz, one needs to decompose the total Hamilto- 
nian of the system into a series of pairwise terms, which 
could be done easily for the model studied by Kurmann et 
al. since all the "pairwise components" of the factorizing 
magnetic field had the same weight in this particular case. 
However, this direct method of determining factorization 
is in fact significantly less obvious, useful, and compre- 
hensive, when one tries to apply it to more general cases. 
For instance, considering models with interaction terms 
of growing spatial range, and/or increasing the spatial 
dimensionality of the lattice, and/or treating extensions 
to spatially anisotropic or non-translationally invariant 
models, results in increasingly hard "guesses" concerning 
the partitioning of the external field into pairwise com- 
ponents and in increasingly nontrivial verification steps. 
Depending on the structure of the model and on the form 
of the product state ansatz, the direct method may ei- 
ther fail to detect factorization points at all, or it may 
detect only a subset of all physically realizable factor- 
ized solutions, or finally it may lead to intractable sets 
of conditions. These limitations of the direct method be- 
come evident, as we will see in the following, even in the 
simplest 1-D XYZ model with nearest-neighbor inter- 
actions analyzed by Kurmann and coworkers. In their 
original work, due to the insufficient generality of their 
proposed product ansatz, they failed to identify a wide 
range of situations that allow the existence of factorized 
ground states. As we will see when discussing models 
with competing interactions of different spatial ranges, a 
further serious fault of the direct method, that becomes 
incurable in the presence of frustration, is its incorrect 
assessment of factorized ground states that are in fact 
excited energy cigenstatcs of the system. 

Following the seminal analysis by Kurmann et at, the 
direct method has been applied to the fully connected 
infinite-range Lipkin-Meshkov-Glick model by Dusuel 
and Vidal [2y|, and to the anisotropic Heisenberg chain 
with the same interactions of arbitrary ran ge a long the 
three spatial directions by Hoeger et al. [24I and by 
Rossignoli et al. [2^ . We would like to remark that these 
works, that include the results of Refs. [23,[2l] as special 
cases, provide an interesting characterization of factor- 
ization points in finite (one-dimensional) translationally 
invariant lattice spin models with periodic boundary con- 
ditions as those points where the parity symmetry is bro- 
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ken. mension and/or with the range of interactions. We will 

show that the difficulties inherent to the direct analytic 
approach and to the numerical method are intrinsically 
B. The quantum informatic approach absent in our formalism. 



A general and rigorous approach to the problem of 
ground state separability, completely different from the 
direct method based on the product state ansatz, has 
been introduced by us in Ref. |27| and is fully devel- 
oped in the present work. This analytic method allows 
inter alia to complete the factorization diagram of the 
anisotropic spin- 1/2 models without restrictions on the 
sign of the couplings, and to go significantly beyond the 
limitations of the previously mentioned studies. Our ap- 
proach is based on an additional characterization of prod- 
uct states (including product ground states) that employs 
concepts and tools of quantum information and entangle- 
ment theory. The strong point of the quantum informatic 
approach to ground state factorization lies in its gener- 
ality and its ensuing independence from the types and 
ranges of interactions present in the Hamiltonian, the 
lattice size, and the spatial dimension, which makes it ap- 
plicable in the most general cases and therefore especially 
useful in all the situations that are either intractable or 
incomplete if one resorts to the direct method or to nu- 
merical techniques. 

In fact, in the original formulation of the quantum in- 
formatic method (27[ there was still trace of a residual 
"ansatz" in the initial stage, since the magnetic order- 
ing of the candidate factorized ground state had to be 
in some sense guessed, and imposed a priori. Moreover, 
in the final stage, we were not yet able to derive a com- 
pletely general set of conditions in order to determine in 
each situation whether a candidate factorized state that 
turned out to belong to the spectrum of the Hamiltonian 
was actually a ground state or belonged to the set of ex- 
cited states. This had then to be decided case by case in 
each specific instance. 

In the present work we proceed further and develop the 
complete method to determine and characterize factor- 
ized ground states in quantum cooperative systems. We 
provide the systematics that allows to determine the oc- 
currence of full ground state separability given a generic 
quantum spin Hamiltonian. The magnetic ordering is 
derived as a consequence of the theory and not imposed 
a priori. A complete set of general conditions is estab- 
lished for the verification of ground state factorization, 
and the elegance of the method does not hinder practi- 
cal usefulness in the application to specific models and 
situations. 

We would like to recall that numerical advances on 
the problem of ground state factorization have been re- 
alized by the Monte Carlo study of ground state factor- 
ization in two-dimensional anisotropic nearest-neighbor 
spin models [2]| . Unfortunately, beyond this important 
result the numerical approach appears to have a very re- 
stricted potential for extensions, since the required com- 
putational power increases rapidly with the spatial di- 



III. SINGLE QUBIT UNITARY OPERATIONS 
AND ENTANGLEMENT EXCITATION 
ENERGIES 



In the present section we provide a brief overview of 
some recent results on the geometrical interpretation of 
the entanglement measure known as "tangle" (linear en- 
tropy) |28| and the consequences it entails for the rela- 
tionship between energy and pairwise entanglement [2^. 
Such results embody the premises for our all-analytic 
characterization of ground state separability in quantum 
spin systems. 

Here and in the following we consider generic systems 
made of N spin-1/2 elementary constituents, or, in the 
language of quantum information, N qubits. For such a 
system, we introduce the set of single-qubit unitary oper- 
ations (SQUOs) Uk [11], defined as the unitary transfor- 
mations that act separately as the identity on any spin of 
the system except the fc-th one, on which they act instead 
as Hermitian, unitary, and traceless operators: 



20, 



(1) 



In Eq. ([T|) li stands for the identity operator on the i-th 
spin of the system while the generic Hermitian, unitary, 
and traceless operator Ok can be expressed as a linear 
combination of the standard spin operators (S^) defined 
on the k-th spin: 

Ok = s'm 0k cos ipkSk +sm6kSm(pkSl + cos OkS^ . (2) 

In the previous definition the parameters 9k and ipk take 
values respectively in the ranges (— 7r/2,7r/2] and [0, 27r). 
Physically, Ok corresponds to a rotation in the spin space, 
and the traceless condition imposes the rotation opera- 
tion to be orthogonal to the identity, so basically a com- 
bination of spin flip and phase flip. 



A SQUO transforms every pure state \'^) G 



^2N 



into a new state \"^) = Uk\'^), that in general dif- 
fers from l^*). We may quantify the action of the 
SQUO in terms of the trace distance between the origi- 
nal and the transformed state, defined as D{Uk', \'^)) = 
D{ek;ipk; |*))=yH(*|C4|¥)|2. The distance D varies 
in the interval [0, 1] and vanishes if and only if the 
two states coincide, meaning that the considered SQUO 
leaves the original state unchanged. For a given initial 
state and for any site k, one can determine the Extremal- 
SQUO (E-SQUO) that minimizes such a distance. One 
can prove that the E-SQUO is uniquely determined by 
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the following conditions on the angular variables [28 



arctan 



— arctan 



{SD cos <p + {SD sin ip 



(3) 



where {S^) denotes the expectation value of the spin 
operator 5'" on the unperturbed state l^*) ((S'^) = 
(\E'|S'^|*)). The E-SQUO plays a crucial role in our anal- 
ysis: For any pure state defined on a general system of 
interacting qubits one can prove [1^ that the square of 



the distance associated to the action of the E-SQUO co- 
incides with the hnear entropy Sl{pu) = 1 — rr[p|] of the 
fc-th qubit (where pk denotes the reduced density matrix 
of qubit k). This quantity, also known as the tangle t in 
the literature, is an entanglement monotone for qubit sys - 
tems and satisfies important monogamy constraints [30j . 
Being a monotonic function of the von Neumann entropy 
of entanglement [3l] | , it quantifies the bipartite block en- 
tanglement present in the state \^) between spin k and 
the remaining N — \ constituents, and is thus a proper 
measure of single-site entanglement [s^l- One has the 
following relation [28j : 



J 



D^Uk. m = ru\y,Mm ^ 4Dctp, = 1-4 [{SD^ + {Slf + {S^f 



(4) 



r 



This geometric-operational characterization of the tangle 
allows to develop a quantum informatic approach to the 
problem of characterizing the correlation properties of 
the ground state in many-body quantum systems [1^. 
Let us consider a collection of spin-1/2 systems defined 
on a regular lattice, with Hamiltonian i?, and let \^) 
denote from now on the ground state of H . If the system 
is in the ground state, the application of a SQUO on 
a generic spin of the lattice is necessarily associated to 
an energy transfer that one must provide to the system 
in order to realize the selected SQUO. In other words, a 
SQUO perturbs the system, and this perturbation results 
in an increase of the average energy (obviously, in general 
the transformed state needs not to be an cigcnstatc of the 
Hamiltonian). Such energy deviation from the ground- 
state value can be quantitatively defined as 

AEiUk) = AEi9k,ipk) = {^\ulHUk\^)-{^\H\^>). (5) 

Because \'^) is the ground state of H, AE{Uk) is a non 
negative defined quantity. Considering the special rela- 
tion between the E-SQUO and the tangle we may expect 
that the associated amount of energy difference AE{Uk), 
appropriately named 'Entanglement eXcitation Energy' 
(EXE), will have as well a strong, direct connection with 
the tangle and in general with the single-site entangle- 
ment in the ground state. In fact, starting from the above 
definition of Eq. ^ and taking into account the property 
of the E-SQUO of leaving unchanged a fully disentangled 
state, it has been established in full generality [1^ that 
if the system is invariant under spatial translation and 
the Hamiltonian of the system does not commute with 
any possible SQUO, then the vanishing of the EXE is a 
necessary and sujficient condition to admit a fully factor- 
ized ground state j^") [2§|. Hence, if the system that we 
want to analyze satisfies all the former hypotheses, the 
vanishing of the EXE provides an additional condition, 
besides the vanishing of the tangle t, that one needs to 



include in establishing a general analytic approach to the 
problem of ground state factorization. 



IV. THEORY OF GROUND STATE 
FACTORIZATION: GENERAL ANALYTIC 
METHOD 

A. Preliminaries 

The general method that was first partially developed 
in Ref. [l^l for the determination of existence, location, 
and exact form of factorizcd ground states rests on the 
following two main observations: 

(i) The ground state j^*) of a spin-1/2 Hamiltonian H is 

factorized if and only if the single-site entanglement 
(tangle) r vanishes for all spins k in the lattice; 

(ii) Provided that [H, Uk] ^ for every possible SQUO 

Uk of the form Eq. ([T]), the ground state \'^) of a 
spin-1/2 Hamiltonian H is factorized if and only 
if the EXE AE{Uk) = {'S\UkHUk\'^) - (^li?!*) 
vanishes Vfc. 



In Ref. [27[ , the strategy to the determination and un- 



derstanding of factorization for a given Hamiltonian bore 
on the above two facts and on the requirement to as- 
sume as working point a phase endowed with some kind 
of magnetic order. By imposing the vanishing of the EXE 
and of the tangle, we showed how to determine uniquely 
the form of the completely disentangled state compatible 
with the assumed magnetic order, candidate to be the 
ground state of the system at a particular value of the 
external magnetic field (labeled as "factorizing field"). 
The final steps of the procedure concern the determina- 
tion of exact analytic conditions for the candidate fac- 
torized state to be an eigenstate of the Hamiltonian with 
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lowest energy. In Ref. [27[; we derived a closed set of 
eigenstate conditions, but establishing whether a factor- 
ized eigenstate is actually a ground state had to be left 
to a case by case analysis. 

In the present paper we generalize the method, extend 
its scope, and apply it to different classes of quantum 
spin models. Firstly, we will show that there is no need 
for a priori assumptions on the magnetic order. In fact, 
we will prove that the existence of a factorized ground 
state is a necessary and sufficient condition for the exis- 
tence of quantum phase transitions and ordered phases in 
translationally-invariant quantum spin models with ex- 
change interactions and in the presence of external fields. 
Moreover, the actual structure of the factorized ground 
state determines automatically the kind of ordered phase. 
Further, we derive a complete set of rigorous conditions 
for the candidate disentangled state to be an eigenstate 
and a ground state. In particular, for systems with no 
frustration we prove that if the factorized state is an 
eigenstate, then it is always a ground state of the sys- 
tem. 



B. Quantum spin models on regular lattices 

We illustrate the method in detail by considering its 
application to the general, translationally invariant, ex- 
change Hamiltonian H for spin- 1/2 systems on a D- 
dimensional regular lattices, with spin-spin interactions 
of arbitrary range and arbitrary anisotropic couplings: 



rr ^ \ ^ jr ax qx , jr qV qV 

2 2-^ ^ i L y i L 

id 



(6) 



Here i (and similarly 1) is a Z)-dimensional index vector 
identifying a site in the lattice, 5" (a = x,y,z) stands 
for the spin-1/2 operator on site i, h is external field di- 
rected along the z direction, r = |i — i| is the distance 
between two lattice sites, and is the spin-spin cou- 
pling along the a direction. Translational invariance is 
ensured by the fact that all the couplings depend only on 
the distance r between the spins. This type of Hamilto- 
nian encompasses a large variety of models and spans over 



several universality classes including, among others, the 
Ising, XY, Heisenberg, and XYZ symmetry classes. Be- 
sides their importance in quantum statistical mechanics 
and in the theory of quantum critical phenomena, these 
models play an important role in the study of various 
schemes of quantum information and quantum commu- 
nication tasks with quantum many-body systems [l^ . 

It is a rather straightforward exercise to verify that the 
Hamiltonian Eq. ^ never commutes with the SQUOs 
Uk Vfc. Translational invariance and non-commutativity 
of the Hamiltonian with the SQUOs guarantee that the 
vanishing of the EXE is a necessary and sufficient con- 
dition for the full separability of the ground state [state- 
ment (ii) above]. Since the external field h is uni- 
form on the lattice, it forces the existence of a site- 
independent, nonvanishing magnetization Mz along the 



M^yk. On the other hand, at a factor- 



z axis: ^ _ 
ization point, statement (i) above imposes the vanishing 
of the single-site entanglement on all sites of the lattice. 
Therefore, if a factorization point exists, the magnetiza- 
tions along the directions orthogonal to the external field 
must assume the following form: 



Mi 



Mj_ cos ipk, 



Ml ^ Mj_ sin (fik, 



(7) 



where M± = y i — A/| is the modulus of the projec- 
tion of the magnetization on the xy plane and the local 
orientations ipk, and hence the type of magnetic order- 
ing in the xy plane, remain undetermined. Since Mz is 
site- independent, then M± is site-independent as well. 
Taking into account Eq. ^ of Sec. IIIII one sees that, at 
least at the factorization point, 9k, i.e. one of the two 
angles that fix the orientation of the E-SQUO, does not 
depend on the site index. Therefore 0^ = 9 Vfc. 

Collecting the former results and considering that in 
the presence of a factorized ground state, if it exists, all 
the correlation functions are products of single-site ex- 
pectations values, we obtain that the condition on the 
vanishing of the EXE associated to the generic site k 
reads: 



J 



AS(c7fe) 



= - sin 61 (sin 9Mz - cos eM±_) Mz ^ J|,, - hf 



■ cos 9M±_ (sin 9Mz ~ cos 9Mi_) I cos ip^ ^ ' Jk.i cos ipi + sin / Jk i 



(8) 



In the former equation /i/ is the factorizing field, i.e. the value, if it exists, of the external field for which groimd 
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state factorization occurs. 

Since we are looking for a fully factorized ground state, 
Eq. ([8|) must be satisfied simultaneously for all k and, 
henceforth, Eq. ^ must be site- independent. This fact 
in turn implies that the following relation must be satis- 
fied by every angle ipk- 

cos ipk Jx cos (fij + sin (pk Jy sin ipj — K . (9) 
L i 

Eq. ([5]) must hold independently of the site on which the 
E-SQUO acts and the value of the constant K must be 
identified by determining the expression of the energy 
density, i.e. the energy per site, associated to the factor- 
ized ground state. Given the total ground-state energy 
E and the total number of sites N , the latter quantity 
can be written as 

r 

+ - M\ ^ COS ipk COS ^pj -I- sin ipk sin ipj 

I 

^ \MlY.^rJ:-hsAU + \MlK, (10) 

r 

where is the coordination number, that is the num- 
ber of sites at distance r from a given site. In order to 
determine the magnetic order in the candidate factor- 
ized ground state, one can now exploit the fact that, by 
definition, it must minimize the energy of the system. 
Therefore such state will be simply characterized by the 
magnetic order that minimizes the per-site contribution 
depending on the Hamiltonian parameters and the 
geometric structure of the lattice. Thus the general con- 
ditions for full separability of the ground state lead to the 
uncovering of the magnetic order in the xy plane. In the 
orig inal formulation of the quantum-informatic method 
[27| one had to assume from the beginning the type of 
ordered phase, and this resulted in constraints on the 
magnetization and limitations in the range of the cou- 
pling coefficients. 

C. Magnetic orders and factorization 

The outcome of the minimization of the energy density 
depends on the presence or the absence of frustration in 
the system, either due to the lattice geometry, or to the 
structure of the Hamiltonian, or to both causes. For 
Hamiltonians belonging to the class of Eq. ([S]), frustra- 
tion can be due both to the geometry of the lattice, e.g. 
for systems with nearest-neighbor antifcrromagnetic in- 
teractions defined on a lattice that admits closed loops 
with an odd number of spins, and to the competition of 
interactions of different spatial orders, e.g. in systems 
defined on linear chains with nearest-neighbor ferromag- 
netic and next-nearest-neighbor antifcrromagnetic inter- 
actions. Frustration effects in systems with interactions 



of different spatial ranges arise whenever it happens that 
the minimization of the energy associated to interactions 
over a certain spatial scale precludes the possibility of 
minimizing the energy associated to interactions over a 
different spatial scale. In the absence of frustration, i.e. 
when the lattice geometry and the interactions between 
the spins are such to allow all energy minimizations si- 
multaneously, it is straightforward to prove that every 
possible factorized ground state must be characterized by 
one of only four magnetic orders out of the many possible 
ones. The complete proof of this statement is reported in 
Appendix IIXI This restricted set includes ferromagnetic 
ordering along the x direction [ipk = OVfc) or along the 
y direction {ipk_ = 7r/2Vfc), and antiferromagnetic order- 
ing along the x direction ((^^ = |fc|7rVA;) or along the y 
direction (ipk = (7r/2) 4- |fc|7rVA:). By requiring the min- 
imization of the energy per site, Eq. pO)) . the type of 
order effectively present in the ground state can be de- 
termined by comparing the values of the different "net 
interactions" along the x and y directions, where the net 
interactions are defined as follows: 

oo 

r=l 
oo 

jf = Y.^rJL, (11) 

r=l 

with a — x,y. The net interaction along the z direction 
is defined equivalently as 

oo 

J'/^^^J^Z.J:. (12) 

r=l 

As a function of the net interactions, the magnetic or- 
ders in the system are determined by the values of 
M = min{jZf,:4^,J-/,J-/}: 

{JjT =4> Ferrom. order along x ; 

Jy Ferrom. order along y ; , . 

J.^ Antiferrom. order along x ; ^ ' 

Jy Antiferrom. order along y . 

It is interesting to observe that, in terms of the net in- 
teractions, one can immediately establish the presence 
or absence of frustration in the system: Any frustration 
arising from the presence of competing interactions would 
imply that ^ — X^r- ^r\Ja\, with a = x,y depending on 
which of the two axes is characterized by a non vanish- 
ing value of the magnetization. We remark that here 
and in the following we are not considering the partic- 
ular situation in which a saturation occurs rather than 
a true ground state factorization: This instance will be 
discussed separately in SubScc. IIVDI 

From now on, without loss of generality, we will spe- 
cialize our analysis to the the case of an antiferromagnetic 
ordering along the x direction. With trivial modifica- 
tions, all the steps and results that will be obtained in 
the following hold as well in the remaining three cases of 
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Eq. (|13p . At the end of the procedure we will present a 
summary table collecting the main results for all the four 
types of magnetic orders compatible with factorization. 

Given an antiferromagnetic ordering along the x direc- 
tion, and keeping in mind that in such case ipk = |fc|7rVfc, 
Eq. d?]) implies that My = while = ±M^ = ±A4 
with = M± > 0. The ± sign reflects the fact that the 
antiferromagnetic order discriminates between two sub- 
lattices, each characterized by the opposite sign of the 
staggered magnetization along x. Consequently, the gen- 
eral condition Eq. ^ takes the form 

(smOM:, ~ cosOMa;) X 

X [cos 9M^J^ - sin 6 (M.jf ~hf)] , (14) 

that admits two solutions for 9, 



tanO = ; tanO — 

JfM,-hf 



(15) 



However, since the E-SQUO must be unique [28[, the two 
solutions must coincide. Factorization thus requires that 



hf 



nrF _ 'jA 



(16) 



Moreover, the requirement of a vanishing tangle (single- 
site entanglement) imposes 



1 



1 - 



Ah) 



iJF-Ji 



A\2 



(17) 



Combining Eqs. p^HTT)) via Eq. ^ yields a closed ex- 
pression for the phase as a function of the Hamiltonian 
parameters and of the factorizing field hf: 



cos 9 



2hf 



^ Z '-' X 



(18) 



Therefore, Eq. ([T5)) . together with the knowledge of 
(fk = T^lfcl, once inserted into Eq. ([2|), determines, in- 
dependently of the value of the magnetizations, the form 
of the E-SQUO at factorization: 



Ok ^ cos 9s:- 



sm9Sf 



(19) 



The form of the candidate factorized state is then 
readily determined by imposing that the action of the 
E-SQUO on all lattice sites leaves the ground state in- 
variant: 



I20fe|*/) 



(20) 



^From Eas. (pni - [TO1) we find then that the candidate 
factorized ground state |^/) has to be a tensor product 
of single-site states that are the exact eigenstates of the 
operators Ok with eigenvalue 1/2: 



I* 



f) 



(21) 



iVfc) = cos(0/2)| Tfc) +e*'^^sin(0/2)| U) . 



Equation ([21]) expresses the general form that a factor- 
ized ground state must assume given an antiferromag- 
netic order along x. This result is independent of the 
actual values assumed by the staggered magnetizations 
Ma's. We remark that in order to establish the gen- 
eral form of the candidate ground state it is essential 
to impose the vanishing of the EXE. In fact, without 
requiring it, the phase 9 and the candidate factorized 
ground state become explicitly dependent on the magne- 
tizations. Then, imposing the conditions for Eq. ([2T|) to 
be an eigenstate of the Hamiltonian, one merely obtains 
a relation between the factorizing field hf and the mag- 
netizations Ma's. If the analytic expression of at least 
one of the Ma's is known, it is still possible to obtain 
an expression of the field hf in terms of the Hamiltonian 
parameters. However, this actually occurs in few special, 
exactly solvable, cases (for instance, the XY model [3^). 
Otherwise, in general one needs to resort to numerical 
evaluations of the magnetizations or introduce other ap- 
proximations in the analysis. This problem is completely 
eliminated in the general analytic framework based on 
the invariance of the factorized ground states under the 
action of E-SQUOs and on the vanishing of the EXE at 
factorization. 

Before proceeding further let us clarify a relevant point. 
Because of the even-odd symmetry of the energy spectra 
of the system, as elucidated also in Refs. [22, 1231, the 
assignment of one of the two different staggered magne- 
tizations to each sublattice is arbitrary. According to 
this observation, one sees that, inverting the sign of the 
magnetization along the a; direction on each site, all the 
conditions for factorization are still satisfied but the can- 
didate ground state takes the form 



(22) 



m = cos(0/2)| tfc) +e*'^i±isin(0/2)| U) . 



Obviously the two states in Eq. (PT|) and in Eq. ([^^ are 
distinguishable and both j^"/) and I'fj) are legitimate 
candidate ground states of the system, as well as gen- 
eral linear combinations of the two (the latter will be in 
general highly entangled). Here and in the following we 
assume to be working in a situation of broken symmetry, 
i.e. after that a small perturbing external field along the 
X direction has lifted the degeneracy between j^*/) and 
I'^'f), e.g. lowering the energy of the former and rising 
that of the latter. After a time long enough to ensure 
convergence to equilibrium, i.e. relaxation to the state 
l^"/), the perturbation is switched off: In the thermo- 
dynamic limit, this ensures that the system will remain 
indefinitely in the state l^*/). 

Having determined the exact form of the candidate fac- 
torized ground state, the subsequent step concerns the 
determination of the conditions for its occurrence, i.e. 
the conditions under which a state of the form Eq. (PT|) 
is indeed the eigenstate of H, Eq. with the lowest 
energy eigenvalue. To this aim, taking the value of the 
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external field at the factorization point, h = hf, it is use- 
ful to decompose the total Hamiltonian as a sum of pair- 
wise Hamiltonian terms Hij, so that H\h=hf = ^jj^ 
where each term Hi j involves only the degrees of freedom 
of a single pair of spins: 

TT jr QX QX I jr qV qV , rr qz qz f,'" / CZ i cz \ 

^hi ^ •Jx^i^j + •-'y^i^j + •■'z^i^j ^ + "-"i J • 

"(23) 

Here the quantity /i^ plays the role of the "component" 
of the factorizing field that acts on the selected pair of 
spins and obeys the following relation: 

2h}^cose{Jl-{-lYJ:) . (24) 

The above ensures that hf = J^,- Zrh^f satisfies Eq. (fTS)). 
Proving that the candidate factorized ground state j^*/) 
is an eigenstate of the total Hamiltonian H is equivalent 
to prove that it is a simultaneous eigenstate of all pair 
Hamiltonians iJy , or, more precisely, that the projection 
of 1^"/) onto the subspace of two given spins i and j — 
which is still a pure state since l^*/) is a tensor product 
of single-spin pure states — is an eigenstate of Hij for 



every pair To proceed, we associate to each lattice 

site a set of orthogonal spin operators defined as follows: 

Ai = cos 6 cos ipiSf - sin 05f ; 
Af = cos ^,Sf- (25) 
= sin 6 cos ipiSi + cos OS- . 

It is immediate to observe that A^ = 0|. By inverting 
Eqs. (pS)) . we can conveniently re-express the standard 
spin operators as functions of the new set of operators 
{Af}: 

S'f = cos ifi (^A^ cos 9 + A- sin 0^ ; 

Sf = Afcosv?,; (26) 
SI = -Af sin6l-t- A|cos6l. 

Inserting Eq. in Eq. we obtain the expression 
of the pair Hamiltonian Hi j as a function of the sets of 
operators {Af} and {A"}: 



H,j = AlA'j [cos^ 6'jr + sin^ 6lj;'(-l)1 - 2/i/ sin 6* 



Af A 



+ A^A'] (sin^ 9Jl + cos^ 0(-l)V:) + A'^Ap;^{-iy - h} cosO (^A^ + A^ , 



(27) 



where we recall that r is the distance between sites i 
and J and that, in the presence of an antiferromagnetic 
order along the x direction, cos ipi cos ipj = (—1)''. The 
spin-pair projection \tpi)\tpj) of IvP/) is an eigenstate of 
the pair Hamiltonian Hij either if all terms in Eq. ([77|) 
admit \^|:i)\■^pj) as an eigenstate or if they annihilate it. 
By construction, recalling that A^ = Oi and that Itpi) 
is the eigenstate of Af with eigenvalue equal to 1/2, the 
spin-pair projection of l^"/) is an eigenstate of the pair 
Hamiltonian if the following condition holds: 



Jy 







(28) 



For general models of the type Eq. ^ with a given, ar 
bitrary. maximum spatial range of interaction s, Eq. ((28 
admits the following solution: 



cos^ e ■ 



(29) 



for all distances r < s, i.e. for all distances associated 
with nonvanishing couplings (we recall that = Va 
and Vr > s -I- 1). The set of equations ([29]) determines 
the conditions that must be satisfied simultaneously by 
the phase 9, which is unique, in order for l^*/) to be an 



eigenstate of the total Hamiltonian Eq. Eqs. ([^^ dis- 
criminate quite clearly between short-range models, i.e. 
models with only one nonvanishing coupling at r = s = 1 
(models with nearest-neighbor interactions) and all other 
possible models containing finite- or infinite-range in- 
teraction terms. Specifically, models with only nearest- 
neighbor interactions are characterized by the fact that 
Eqs. ([29|l reduce to a single condition and hence if the 
value of the r.h.s. of Eq. (^9]) falls in the interval [0, 1] (as- 
sociated to the permitted values of cos^ 9 in terms of the 
interaction parameters), one can immediately conclude 
that the models under investigation admit a factorized 
eigenstate. On the other hand, for models containing 
interaction terms of longer range — so that Eqs. (|29|) in- 
clude two or more conditions — factorized energy eigen- 
statcs are allowed if and only if the r.h.s. of all conditions 
in the set ((29)) take values in the interval [0, 1] and, more- 
over, they all coincide. 

Summing term by term over the index r all the rela- 
tions in Eq. ((^ . taking into account Eq. ((T5() . and solv- 
ing for hf, one eventually obtains the exact expression of 
the factorizing field as a function of the net interactions: 



2\/('^/ ~ ^z){^^ ~ ^z) 



(30) 
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We remark that Eq. pO|) , which was derived in Ref . [27[ 
using some unnecessary auxiliary assumptions, is com- 
pletely general and holds for lattices of arbitrary spa- 
tial dimension and for spin-spin interactions of arbitrary 
range. Accordingly, the angle 6 that determines the di- 
rection of the E-SQUO can be also expressed as a func- 
tion of the net interactions as follows 



(31) 



In the case of a maximum range of interaction s > 2, 
relation Eq. ([?T|) for the net interactions is a necessary 
but not sufficient condition for factorization, and further 
use of Eqs. (|28|) and Eqs. (|29|) is needed, as we will show 
with explicit examples in Sec. IVII and Sec. IVIlIl 

We are finally left with the problem of establishing 
conditions for a factorizcd energy eigenstate j^*/) to be 
indeed a ground state of the system. A very simple suf- 
ficient but not necessary condition is that the spin-pair 
projection of onto the subspace of a pair of spins be 
the simultaneous ground state of every pair Hamiltonian 
Hij. In fact, given the relation existing between all the 
pair Hamiltonians Hjj in Eq. and the total Hamilto- 
nian H Eq. ([6]), it follows that the simultaneous ground 
state of all pair Hamiltonians Hjj is also ground state 
oi H. To proceed, let us re-express the pair Hamiltoni- 
ans Eq. ()27p in matrix form in the basis spanned by the 
eigenstates of the operators (E) Aj. One has 





—a,. 
J^/2 





.rj2 



\ 





A- 

-(-i)'^A 



(32) 



where = ^{cos^ej'^ + (-l)" sin^ ) and Pr = 
h'j: sm0. Making use of Eq. and Eq. ([M]) the eigen- 
values of the matrix Eq. (|32[) can be put in the form 

so = I [3(-l)Vi - Jr - ( JI - (-l)V^) COS (20)] 
ei = i[3(-l)Vi- Jr-(JI-(~l)VDcos(20)] 
e2 ^ I [-3 {{-lyj^. + JI) + ( JI - (-l)VD cos (2^?)] 

es ^ I [-H-iyr, + 5jr + ( j: - {-lyrj cos {29)] 

(33) 

The eigenvalue Eq is the one associated to the projec- 
tion \ipi)\ipj) of the factorized eigenstate j^"/) onto the 
four-dimensional Hilbert space associated to the pair of 
spins {5*^,5^}. Therefore, the spin-pair projection is the 
ground state of the pair Hamiltonian if Eq < {ei, £2, £3}- 
Imposing this condition yields the following inequalities: 



Notably, these inequalities are automatically satisfied for 
any interacting spin system that is not frustrated and 
that verifies all of the Eqs. ([29|) . Namely, for every even 
value of the distance r in the interaction range (r < s) 
Eq. (|35[) implies that the denominator in the correspond- 
ing Eq. ((29|) is non-positive. Therefore, an acceptable 
value of cos'^ 9 can be obtained only if also the numera- 
tor is non-positive and exceeds or equals the denomina- 
tor. These conditions in turn imply the order relation, 
JI < Jr. On the other hand, comparing Eq. ([55)1 and 



ij: + jy)[j:^i~iyj:] > o 
{~iyj:-j: < o 



(34) 
(35) 



Eq. fM]) . one has that < —Jy- Therefore, compar- 
ing the two order relations yields finally < —\Jy \ {r 
even). This condition is always verified in systems with 
an antiferromagnetic order along the x direction and in 
absence of frustration and therefore the projection of the 
factorized energy eigenstate on the Hilbert space of ev- 
ery pair of spins is indeed the ground state of every pair 
Hamiltonian. On the other hand, in the case of odd r, 
Eq. ((55l) implies that the denominator of Eq. is non- 
negative and, similarly to the former case, the numerator 
must be non-negative and not exceeding the denomina- 
tor. Hence we obtain > Jy while from Eq. ([M)) we 
recover > — that, in turn, implies > \ Jy\. Again, 
such inequalities are always verified in the presence of an 
antiferromagnetic order and in the absence of frustration. 



Collecting all these results, we have proved that the 
following holds: 

Theorem. For any cooperative system of spin-l/2 par- 
ticles described by Hamiltonians H of the type [Eq. 0)/, 
characterized by an anti-ferromagnetic order along the x 
axis as emerging from Eq. and in the absence of 

frustration, the simultaneous verification of all Eqs. i2^) 
is necessary and sufficient for the fully factorized state 
Eq. \21\) to be the exact ground state of H when the ex- 
ternal magnetic field takes the value h = hf determined 
by Eq. 

This central result holds as well, with obvious modi- 
fications, when one considers the other different ordered 
phases that form the set identified in Eq. In Tab.|T] 

we summarize and compare results for each of the four 
different possibilities. 



D. Factorization, balancing, and saturation 

Before moving to apply the general method to specific 
models and examples of conceptual and physical rele- 
vance, we discuss briefly the meaning of factorization and 
its marked differences with the phenomenon of "satura- 
tion" (see below) . As we have seen, and as it will appear 
even clearer in the discussion of the examples in the fol- 
lowing sections, the physical mechanism of ground state 
factorization in quantum spin systems is due to a deli- 
cate kind of "balancing" between the coupling strengths 
that regulate the intensity (and range) of the interactions 
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TABLE I: Summary of results for the four magnetic orderings [Eq. (|13|) ] compatible with ground state factorization. 



magnetic order cos^ 9 factorizing field hf single-spin state \ipi^ 



ii^Jl cos(e/2)| Tfe)+sin(e/2)Ufc) 



M = :^/ \^{J^-JFW-JF) cos(e/2)| T.) + (-1)1^1 sin(e/2)| U> 



ii = J^ T^f \^{Ji -JF){JF -3F) cos(e/2)i U) + ^sin(e/2)| U> 



A' = ^/ 'if^^S'h h^i^^^-^FW-JF) cos(e/2)| Tfc> + (-l)l^lisin(e/2)| U) 



and the aligning effect of the external field, that tends to 
orient all spins along a given direction and to destroy all 
quantum correlations. The remarkable aspect of factor- 
ization is then that it occurs, according to well defined 
conditions and constraints, at precisely defined, finite val- 
ues of the couplings and of the fields in an ordered (or 
symmetry-broken) phase. This means that factorization 
occurs when the system is relatively strongly interacting 
and the external aligning field is relatively weak. More- 
over, we will see that the factorization point is always 
a precursor from below (from the ordered phase) of the 
quantum critical point, as first observed for the XYZ 
model with nearest-neighbor interactions [2l[ . 

Saturation is the phenomenon of ground state factor- 
ization that occurs trivially when the value of the exter- 
nal field grows unboundedly compared to all other Hamil- 
tonian parameters. Therefore its aligning effect on the 
spins prevails against all other effects and quantum cor- 
relations are suppressed. The main physical difference 
between proper factorization and saturation is then clear 
and is reflected at the level of the Hilbcrt space of states. 
Namely, suppose that the system under study admits a 
fully factorized state that is the true ground state of the 
system at h = hf. Then, when one moves away from the 
factorizing field, i.e. when h — hf = e =^ 0, the factorized 
state ceases to be an eigenstate of the Hamiltonian. On 
the contrary, in the presence of saturation, the trivially 
ensuing factorized state is always an eigenstate of the 
Hamiltonian and remains the ground state of the system 
when perturbing the value of the external field h. 

Let us formalize the above discussion. Consider a sys- 
tem that admits a fully factorized ground state. In this 
case the pair Hamiltonian in Eq. I|27p admits as eigen- 
state a factorized spin-pair state. However, Eq. P7|) is 
the expression of the pair Hamiltonian only at h ~ hf. 
If we consider a value of the external field different from 



the factorizing field we must include in Eq. (P7|) a term 
proportional to {Af — A^)sm9{h — hf). Obviously, for 
nonvanishing values of sin 9 the presence of such term pre- 
vents the factorized state from being an eigenstate of the 
Hamiltonian of the system. On the contrary, if sin 9 — 
and hence if Ok = Vfc the extra Hamiltonian term van- 
ishes and the factorized state {^j, | 1)^ remains an energy 
eigenstate of the system for all values of the external field. 
This result can be understood from a different point of 
view. When 9 = 0, from Egs. ipS)) we immediately ob- 
tain that = Jy, which implies that the total magne- 
tization along the z axis is preserved, [H, Sf,] = 0. 
Therefore the Hamiltonian of the system and the total 
magnetization operator admit a complete set of common 
eigenstates, discriminated by the eigenvalue of the total 
magnetization: the factorized state is the one (and only 
one) characterized by the maximum eigenvalue. 

We now move on to show that the factorized eigen- 
state, in the case of saturation, is the system's ground 
state in an infinite range of values of the external field. 
Here and in the following we recall that since the exter- 
nal field is directed along the conventional z direction, 
only saturation with spin alignment along the z-axis is 
allowed. At very low values of h the ground state of the 
system will possess a small total magnetization along the 
field direction. However, as h increases, also the total 
ground-state magnetization along the z direction will in- 
crease until all spins align with the field and saturation 
occurs. By increasing h further, it is not possible for the 
system to evolve in a state with higher values of the mag- 
netization along the field direction. Hence, the factorized 
state remains the ground state of the system for all values 
of the external field h > hs, where by hs we denote the 
value of the field for which the system reaches saturation 
("saturating field"). Let us consider the matrix represen- 
tation, similarly to Eq. (|32p . of the pair Hamiltonian at 
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arbitrary values of h. Taking into account the condition 
sin 6 = we have 



/ JJ/4 - h'^ 



V 





Jill 






JII2 




(36) 

It is immediate to verify that if > | J^|/2 + J!72 = h^, 
the projection of the factorized state onto the consid- 
ered pair of spins is associated to the lowest energy. Re- 
summing over all values of the spatial range index r and 
taking into account the coordination number Zr of the 
lattice, we have that the factorized state is the ground 
state of the total system if 



h > h. 



Zr\j:\ 



ZrJl 



(37) 



Eq. p7p is the rigorous condition for saturation. We 
notice that it coincides with Eq. ((30|l in the limit Jy 
J^, meaning that in this simple instance no ground-state 
factorization occurs other than saturation. 

In the following Sections we apply the general theory of 
factorization to different quantum spin models. We will 
derive a series of exact results concerning the occurrence 
of factorized ground states in systems with various types 
of two-body interactions of different spatial ranges de- 
fined on regular lattices of arbitrary spatial dimensional- 
ity. The general method enables to determine many novel 
exact factorized ground state solutions as well as the ex- 
istence and the nature of ordered phases in models that 
are, in general, non exactly solvable. Moreover, it allows 
to recover quite straightforwardly the existing analytical 
and numerical results. In the process, we will come to ap- 
preciate that ground state factorization is a phenomenon 
more common than previously believed. Determining the 
existence of exact factorized solutions in non exactly solv- 
able models should in turn allow to envisage controlled 
approximation schemes, e.g. perturbative or variational, 
to investigate the physics of non exactly solvable quan- 
tum spin models in the vicinity of quantum factorization 
points. 



V. MODELS WITH SHORT-RANGE 
INTERACTIONS 

In this and in the following sections we will consider 
quantum spin models of the XYZ type Eq. ^ defined 
on regular lattices of arbitrary dimension. Let us first 
consider models with short-range interactions including 
only nearest-neighbor couplings: 



dimension, whose geometry does not induce any frustra- 
tion, i.e. in which any closed loop is characterized by an 
even number of spins. The absence of geometrical frus- 
tration together with the fact that all the couplings be- 
yond the nearest neighbor spins vanish, ensure that it is 
always possible for the system to satisfy all the conditions 
on the net magnetic interactions. Therefore, as shown in 
SecHVl if the r.h.s. of Eq. — or, better, if the r.h.s. 
of the equation relative to each possible magnetic order 
(See Tab. IJ) takes an acceptable value (0 < r.h.s. < 1) — 
then the system admits a fully factorized ground state. 
In the case of antiferromagnetic ordering along the x di- 
rection, which is obtained if Ja; > \Jy\, from Eq. ([29)1 it 
follows that a factorized ground state exists ii Jy < —Jz- 
Analogous conditions are obtained for the other types of 
magnetic orders compatible with full ground state sepa- 
rability. The results are summarized in Fig. [l] where Jz 
is set assume positive values for reference. A very similar 
diagram of factorization holds in the case Jz < 0. 

As one can see from Fig. [U the type of magnetic order 
corresponding to the factorized ground state determines 
four different regions in the space of the Hamiltonian 
parameters (one for each different ordering) associated 
to the existence of a factorized ground state. The four 
different domains of the interaction parameters are sepa- 
rated by two bisectrices that individuate special cases. 
The line J^ = Jy is constituted by the set of points 
in which saturation occurs in place of genuine factor- 
ization (see SubSec. |IVD[) . On the other hand, the line 
Jx = — Jy, —^<Jx/Jz<^, corresponds to models that, 
for a given value of the factorizing field hf and of the 
E-SQUO orientation 9, admit two degenerate, but dis- 
tinct factorized states. The choice of the ground state 
is then determined by the usual mechanisms of symme- 
try breaking: If an infinitesimal magnetic field is applied 
in the xy plane, the axis of magnetic alignment in the 
ground state is determined by the order in which the two 
limits hx ^ Q and hy are realized. 

It is important to observe that all the shaded regions of 
Fig. [1] correspond to exact analytical ground-state solu- 
tions of the investigated short-range models, which are, in 
general, non-exactly solvable. The exact, fully factorized 
ground states are tensor products of the form reported 
in the fourth column of Tab. [H in which the value of 
the orientation parameter 9 of the E-SQUO is retrieved 
by solving the equation in the second column of Tab. [D 
Obviously, such an exact solution holds only at the fac- 
torizing field h = hf. In the case of antiferromagnetic 
order either along the x or the y direction, the latter is 
given by (See Tab. IJ) 



(39) 



Ji = j„, J^ = 0,Vr>2. 



(38) 



while in the case of a ferromagnetic phase, again unre- 
spcctfuUy if along x or along y, it reads 



As in Sec. IIVI we restrict ourselves to the case in which 
such models are defined on regular lattices of a generic 



Z I 

hf ^ ^\l{Jx- Jz) {Jy - Jz) ■ 



(40) 
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FIG. 1: (Color online) Diagram of ground state factorization 
as a function of the ratios Jx/Jz and Jy/Jz- The coupling 
Jz is set to be positive for reference. Analogous results hold 
in the case Jz < 0. Blank regions: No factorization allowed. 
Rightmost (yellow online) shaded region: Presence of a fully 
factorized ground state supporting an antiferromagnetic or- 
der along the x direction. Leftmost (light blue online) shaded 
region: Fully factorized ground state with a ferromagnetic 
order along x. Bottommost (light red online) shaded region: 
Fully factorized ground state with a ferromagnetic order along 
the y direction. Topmost (light green online) shaded region: 
Fully factorized ground state with an antiferromagnetic order 
along y. The solid line — Jy corresponds to models that 
exhibit saturation rather than true factorization; the dashed 
line Jy = —Jx accommodates for models with unresolved su- 
perpositions of two fully factorized ground states each one 
supporting a different magnetic order. Either one of the two 
possible ground states is then "chosen" by the system depend- 
ing on the way the symmetry is broken, see text for further 
details. All plotted quantities are dimensionless. 

Given that the ground state is known exactly at the fac- 
torization point, all the relevant physical observables can 
be evaluated straightforwardly. For instance, the ground 
state energy density at factorization reads 

§ = -Y^'^- + '^y + '^^^ (41) 

in the presence of an antiferromagnetic ordering, while 
in the ferromagnetic case it is 

I = Y ^^^^ + " -^'^ • ^^^^ 

Both the factorizing field and the ground state energy, 
unlike the E-SQUO orientation 9 and the expression of 
the factorized ground states, are functions of the num- 
ber of nearest neighbors per site (the coordination num- 
ber Zi, where the index 1 denotes the nearest-neighbor 



range of the interactions). Hence, they depend on the 
geometry and the spatial dimension of the lattice. To 
begin with, let us consider one-dimensional lattices (spin 
chains, D = 1), for which the coordination number 
Zi = 2. Substituting this value in the previous equa- 
tions, we immediately recover, among others, the value of 
the factorizing field and the ground state energy obtained 
analytically by Kurmann et al. for the XYZ spin chain 
using the direct method [1^. Kurmann et al. had to 
restrict their investigation to coupling constants Jq, tak- 
ing values in the interval [0, 1], thus establishing ground 
state factorization only in a restricted set of conditions. 
Instead, the general method allows to assess the occur- 
rence of factorized ground states associated to different 
types of magnetic orderings in the entire region of the 
Hamiltonian parameters. 



One of the most important advantages of the general 
method is that it is easily applicable irrespective of the 
spatial dimensionality, while the direct method becomes 
increasingly difficult to apply with growing spatial di- 
mensions of the lattice. Numerical techniques can over- 
come some of the limits of the direct method, although 
they fall short of the power and extension of the gen- 
eral analytic method. The most important numerical 
result has been obtained by applying quantum Monte 
Carlo techniques to the study of the XYZ model in 
D = 2 space dimensions (square lattice) [2l|. Apply- 
ing the present method to the XYZ model on the square 
lattice allows not only to immediately reobtain analyti- 
cally the previous numerical findings, but also, again, to 
extend the investigation to the entire space of the Hamil- 
tonian parameters by simply taking into account that for 
the square lattice Zi = 4. Moreover, the flexibility of the 
general method allows to include in the analysis different 
planar models. For instance, considering models defined 
on planar hexagonal lattices, the previous analysis needs 
to be modified only in that the value of the coordination 
number, in this case is Zi = 3. 



Let us now move to D = 3 and investigate ground state 
factorization in nearest-neighbor XYZ models on three- 
dimensional regular lattices, for which there were so far 
no results (apart the preliminary ones obtained in our 
previous work [2^ ) either analytic or numerical (such a 
task would be extremely demanding if approached either 
via the direct method or via numerical techniques). Re- 
sorting to the general method, we can instead fully char- 
acterize the factorized ground state and the associated 
physical quantities just by entering the correct value of 
the coordination number, e.g. Zi = 6 for a simple cubic 
lattice. Obviously, the results extend straightforwardly 
also for models defined on lattices of arbitrary higher di- 
mension D > 3. 
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VI. MODELS WITH FINITE- AND 
LONG-RANGE INTERACTIONS 

The analytic method can be apphed with minor com- 
phcations to models with interactions of longer spatial 
range. As already mentioned in Sec. IIVI concerning 
ground state factorization, the main difference between 
models with short- and finite-range interactions is that 
for the latter there exist more than one equation in the 
set p9)) . Each of these equations is associated to a dif- 
ferent value of the range r < s, where s is the maximum 
interaction range, that is the maximum distance between 
two directly interacting spins: All the couplings vanish 
for pairs of spins with inter-spin distance greater or equal 
to s + 1. 

Clearly, the fact that all the equations have to be solved 
simultaneously yields tighter conditions on factorizability 
compared to the case of models with nearest-neighbor 
interactions. Namely, it is not sufficient that the r.h.s. of 
each Eq. ([M)) takes values in the interval [0, 1]. It is also 
necessary to require that all these values coincide. This 
constraint yields a set of further conditions that must be 
satisfied by all the nonvanishing couplings. For instance, 
let us suppose that the system is in an antiferromagnetic 
ordered phase along the x direction. Labeling by c the 
value of cos^ 9 and imposing that for every r the r.h.s. of 
Eqs. ((29|) must equate c, we find that 

j;; = cr, + {~ini-c)ji w<s, (43) 

where, to ensure the absence of frustration, must obey 
the following inequality: (-l)V^ < JJ(c-l)/(l+c) Vr < 
s. Analogously, for an antiferromagnetic order along y: 

J:^cJl + {-ini-c)j: yr<s, (44) 

and to exclude frustration we must also have that 
i-^yjy < Jl{c + c) Vr < s. In the presence of a 

ferromagnetic phase, respectively for an ordering along x 
and along y, the constraints and the no-frustration con- 
dition become 

JI = cJ: + {1-c)J: J:<'y^J: Vr<s;(45) 

= cJ,'; + (1 - c) jr<'--ljl Vr<5.(46) 

Eqs . generalize the instance analyzed in Ref. 

[27j . in which we limited the analysis to the case of 
vanishing that, according to Eqs. (|43H46|) . implies 
c= J^/J^ for all r < s. 

Especially in the presence of a ferromagnetic order in 
the lattice, Eas. (l45l - [46|) allow us to simplify the expres- 
sions of both the factorizing field and the energy per site. 
Namely, taking into account the definition of the ferro- 
magnetic net interactions Eq. (|11[) and the expression of 
the factorizing field in the presence of a ferromagnetic 
order (see Tab. H]), wc have 

hf = V~c{Jy^-j!) , (47) 



respectively for an order along the x and the y directions. 
Similarly, for the energy density we obtain 

| = (l + c)^/-cj'/, 

|-(l + c)J'/-cJ'/. (48) 

We now discuss a specific example (out of many pos- 
sible ones) which demonstrates the power of the analytic 
method in conditions where the more conventional direct 
approach based on the factorizcd ansatz of Refs. [13, 
is completely ineffective without the supplement of a re- 
markable intuition in guessing the correct forms of the 
possible factorized ground states. Let us consider a 
XFZ-type spin model defined on a spatially isotropic cu- 
bic lattice of arbitrary size, possessing nearest-neighbor 
and next-nearest-neighbor interactions, and admitting 
arbitrarily different couplings along the x, y, and z direc- 
tions. The coordination numbers associated to the two 
types of interaction can be straightforwardly computed 
from the geometry of the problem and are Zi = 6 and 
Z2 = 18. Suppose the model is deduced from an ex- 
plicit experimental realization of a spin lattice and we 
are provided with explicit values of the various interac- 
tion strengths. For example, = 1' "^y = 0.3, = 
0.4, J2 ^ -0.6, J2 = -0.25, J2 = 0.1. The application 
of our "on-demand" ground-state factorization search en- 
gine is now immediate. From Eq. (|13p we have that 
the candidate factorized state must possess antiferromag- 
netic order along the x direction. Therefore the general 
conditions for the existence of a factorization point, as 
determined by Eq. ((29)) . can be summarized as follows 

cos^ Q — '^y ^ '^^ _ "^y "^^ _ "^y ^ '^z ^ 3( ~ Jz) 
Jx + "^i '^y^ '-^z Jx ^ Jz ^ 3( — ) 

(49) 

This condition must hold together with the frustration- 
free constraint. It is readily verified that cos^ — Q.h pro- 
vides a solution to the combined set of necessary and suf- 
ficient conditions for ground state factorization. There- 
fore if the experimentalist tunes the external field h at the 
factorizing value hf =~ 7.425, obtained from Eq. (|30p . 
the engineered cubic spin lattice system will relax in a 
completely factorized ground state of the form given in 
Table U (second row, last column) with 9 = arccos(-\/0.5). 

Many other applications can be considered to models 
with arbitrary spin-spin interactions of longer range, and 
they can be treated with the same simplicity within the 
framework of the general analytic method. For instance, 
a very interesting example of a model with infinite-range 
interactions will be treated in the next Section. However, 
before moving to the study of this case, we should com- 
ment on the important issue of the interplay between fac- 
torization and frustration in quantum spin systems with 
multiple spatial scales of interaction. 

The investigation of general models including many in- 
teraction terms of different spatial ranges needs to be car- 
ried out with particular care, because the interplay be- 
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tween the different interactions can lead in general to im- 
portant effects of frustration that tend to suppress the oc- 
currence of ground state factorization. The very delicate 
and tricky nature of the problem is at the origin of a re- 
cent incorrect prediction about the existence of factorized 
ground states in one-dimcnsional XY Z models on finite 
rings with long-range interactions, obtained exploiting 
the direct method jsj. In Ref. [s^j the author, resorting 
to the direct method, concludes that all the factorized 
energy eigenstates in one-dimensional long-range XY Z 
models, with interaction terms of spatial ranges that are 
in integer ratio with the total number of sites, are also the 
states of lowest energy. However, the proof of this state- 
ment obtained using only the direct method requires that 
there must be no frustration in the system (20j . Unfortu- 
nately, the interplay between the interactions of different 
spatial range can in general introduce frustration effects 
(e.g. in the case of antiferromagnetic nearest-neighbor 
and next-to-nearest neighbor interactions, and so on). 
The existence of these effects nullifies all proofs, based 
on the direct method, of the existence of fully separable 
ground states. In fact, immediate counterexamples can 
be given of fully factorized energy eigenstates that are 
not ground states. Let us consider for instance an XY 
model on a finite ring of six spins with nearest-neighbor 
and next-to-nearest-neighbor antiferromagnetic interac- 
tions, and Jl = 1, Jl = 0.3, Jl = 0.6, = 0.18. Fol- 
lowing Ref. [3J|, this model should admit a factorized 
groimd state. In fact, a factorized energy eigenstate with 
antiferromagnetic order does indeed exist, with an energy 
E = —0.78. However, solving the model by exact diag- 
onalization yields that the lowest energy is E ~ —1.11 
and that the associated ground state is entangled. Es- 
sentially the same type of counterarguments nullifies a 
recent claim that the factorized energy eigenstates are 
the factorized ground states of mixed-spin models of fer- 
rimagnetism of the XYZ type [IHl ■ A general analysis of 
ground state factorization in ferrimagnetic models with 
mixed smn- 1/2- spin- 1 interactions will be presented else- 
where [36| , based on a proper extension of the formalism 
of SQUOs to include spin-1 systems, as sketched in the 
final part of Ref. 

In conclusion, the study of factorization in quantum 
spin models with many different finite- and long-range 
interaction terms requires a very careful analysis and the 
use of the general analytic method whenever frustration 
effects are present. A rigorous and systematic study of 
the crucial (and subtle) interplay between frustration and 
factorization in quantum spin systems will be the subject 
of a forthcoming paper [33| • 



VII. LIPKIN-MESHKOV-GLICK MODEL 

A very interesting limit of the cases analyzed in Sub- 
Sec. IVII is given by models with infinite range interac- 
tions, such as the fully connected or Lipkin-Meshkov- 
Glick (LMG) model [33 , extensively studied in con- 



densed matter physics [2^, [H, H^l, which is described 
by the following Hamiltonian 

z,j i 

(50) 

with < A < 1. The pre-factor l/4(iV - 1) (with N 
denoting the total number of spins) ensures the linear di- 
vergence of the ground state energy as N increases. Let 
us remark that normalizing the interaction part of the 
Lipkin-Meshkov-Glick Hamiltonian requires some care: 
Being a connected interaction, it scales as x [N — 1), 
i.e. as the product of the total number of spins and the 
number of pairs of spins. Therefore, to ensure that the 
interaction energy is extensive, i.e. scales with N , the 
interaction part of the Hamiltonian must be normalized 
by N — 1. Usually, for calculational ease, the interac- 
tion part is normalized by A^, so that it scales as — 1. 
Using this convention is strictly speaking incorrect, but 
is clearly harmless for large N and in the thermody- 
namic limit. For instance, concerning the separability 
of the ground state, it introduces a slight, spurious A^- 
dependence of the factorization point through the multi- 
plicative factor {N — 1)/N, and the error vanishes in the 
thermodynamic limit. However, on fundamental phys- 
ical grounds, the factorization points of translationally 
invariant systems, if they exist, must be independent of 
the size of the system, and must take the same value both 
at finite A^ and in the thermodynamic limit. Applying 
the Kurmann-Thomas-Miiller factorized ansatz [20| and 
the 1/A^ normalization, Dusuel and Vidal obtained an 
analytic expression of the factorizing field for the LMG 
model that is A^-dependent through the multiplicative 
factor (A^— 1)/A^ [2y|. Although defined according to an 
anomalous convention, this result is essentially exact in 
the limit of large A^. 

We first observe that the LMG model in Eq. (|50p can 
be obtained as the limiting case of a generic model de- 
scribed by Eq. ^ assuming that the range of interactions 
diverges (i.e s — > oo) while the individual strengths be- 
come independent of the distance (J^ = Ja), scale with 
A^ (Jq, (X 1/A^), and vanish along the z-direction ( = 0). 
Therefore, assuming J^' = —2/N, Jy = —2A/N and 
J2 = in Eq. ^ we recover Eq. ([501) . 

Because of the sign of the interactions and of the range 
of values that can be assumed by A, we immediately 
obtain that the energy per site IjlOp is minimized by a 
state possessing ferromagnetic order along the x direc- 
tion. Moreover, taking into account that the coupling 
strengths are proportional to the ratio 1/A^, we have 
that the net interactions converge to finite values in the 
thermodynamic limit: J^f —2, J'^ — > — 2A. It is 
then rather straightforward to apply the analytic method 
to the LMG model and hence prove, exactly, that the 
ground state is fully factorized at the factorizing field 
hf — a/a, which is, correctly, independent of A^. The 
form of the fully separable ground state is the tensor 
product of local states of the form \tpk) = cos(6'/2)| '\k 



) + sm{9/2)\ Ik) with 6 = arccos(-\/A) and the energy 
density reads E/N = —(1 + A)/4. These results are in 
agreement with those obtained, using the direct method, 
by Dusuel and Vidal [2^ , modulo a proper normalization 
factor. 



VIII. MODELS WITH SPATIAL 
ANISOTROPIES 
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The versatility and generality of the analytic method 
can be demonstrated further by considering extensions 
of the translationally-invariant models discussed in the 
previous Sections. In fact, till now, we have always con- 
sidered systems in which the coupling amplitudes depend 
only on the distance between the spins involved. In the 
present subsection we shall focus instead on models which 
contemplate the possibility of spatial anisotropics, i.e. 
given the spin 5*^ associated to site i, the interaction cou- 
pling with spin 5"^ associated to site j depends not only 
on the distance r — \i— j\ but also on the location of site 
j relative to site i. 

To simplify the notations, the physical analysis, and 
the number of possible cases that need to be consid- 
ered, we limit ourselves to models with short-range in- 
teractions, i.e. with all the couplings vanishing for 
r = \i — j\ > 2. We will comment briefly on general 
models with spatial anisotropies and interactions of ar- 
bitrary range at the end of this Section. Let us assume 
that for each site i the total number of nearest neigh- 
bors is Zi = z[^^ + z'i^\ where Z^^'' is the number of 
nearest neighbors whose coupling strength with Si takes 

(2) . ~ 

a certain value fa, and Z^ is the number of nearest 
neighbors whose coupling with Si is given by a different 
value ga- Obviously, we could consider situations more 
complicated at will, with an arbitrary number of n dif- 
ferent types of nearest neighbors. For the given example, 
the system Hamiltonian reads 



H = 1/2^ 



]f.s^si +fySfSl +hS^S, 



^kl 



J2 9.SISI + gySfSl + g.StSl - 2hSl 



(51) 



where the index runs over the first z'^^ nearest neigh 

bors, and the index runs over the remaining Z^ 
nearest neighbors for each spin Si^. According to the 

different values of z\^^ and this Hamiltonian de- 
scribes various possible models. For instance, if Z^^' = 2 

and zf'^ = 1 Eq. ((5T|) describes XYZ models on lad- 
ders constituted by two coupled linear chains; the case 
z'"^^ — 4 and zf^ — 1 corresponds to models defined on 
two coupled square lattices. In both cases, every spin in 



(2) 



FIG. 2: (Color online) Summary of the possible magnetic ar- 
rangements along a single direction [x or y) characterizing the 
candidate factorized ground states of models with spatially 
anisotropic interactions. The case considered corresponds to 
a ladder of two coupled one-dimensional spin chains with co- 
ordination numbers = 2 and Zj^' = 1. Topmost row: 
Parallel (a) and antiparallel (b) ferromagnetic order. Bottom- 
most row: Parallel (c) and antiparallel (d) antiferromagnetic 
order. 



each chain (plane), besides the usual interactions within 
a single chain (plane) involving z'^^ nearest neighbors. 



.(2) 



companion spins located on the 



The instance zf ' = zf ^ = 1 



is also coupled to Z| 

second chain (plane) 
corresponds to a XYZ-type model defined on a linear 
chain with nearest-neighbor interactions of alternating 
strengths. The factorization properties of this latter 
model have been studied by G. L. Giorgi in the limit 
/z = 5z = that, belonging to the XY symmetry class, 
can be solved exactly 13J| . 

As we have seen in Sec. lIVi the first step of the method 
is to single out the magnetic order that minimizes the en- 
ergy associated to the candidate factorized ground state, 
at fixed magnetization along the z direction (the direc- 
tion of the external field). Provided there are no effects 
of frustration, in the present case we obtain that every 
possible hypothetical factorized ground state can assume 
one of eight different magnetic orders. These eight pos- 
sibilities stem from the four possible orders of Eq. (|13p . 
that the system can support when the effects of ga can 
be neglected, duplicated according to the two possible 
arrangements, parallel or antiparallel, induced by the in- 
teractions associated to ga- The four possible orderings 
along a single direction are represented in Fig. [5] for a 
model of two transversally coupled one-dimensional spin 
chains. 

By imposing the minimization of the energy per site, 
Eq. (jlOp , the kind of order actually present in the ground 
state can be determined comparing the antiferromagnetic 
and ferromagnetic net interactions along the x and y 
axes, defined as 



.(1) 



^1 ' fa^ Z[ ^ ga 
.(1) 



Ja ~ +-^1 fa ^ Z[ ^ ga , 

where the ± sign discriminates between paral- 
Icl/antiparallcl ordering with respect to ga- As 
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a function of the net interactions, the type of 
magnetic order present in the system is deter- 
mined by the value of the minimum, i.e. of /i = 
mini 7^^ 7^" 7^* 7^" 7^"" 7^" 7^^ 7^" I 
Accordingly, we have the following correspondences: 





4> Parallel ferromagnetic order along a; 




^ Antiparallel ferrom. order along x ; 


^C- = 


^ Parallel ferrom. order along y ; 


Jy ^ - 


^ Antip. ferrom. order along y ; 


Jx ~ 


Parallel antiferrom. order along x ; 


Jx ^ ~ 


4> Antip. antiferrom. order along x ; 


Jy 


4> Parallel antiferrom. order along y ; 




4> Antip. antiferrom. order along y . 



(52) 



Obviously, also in the presence of spatial anisotropics we 
may observe that, in terms of the net interactions, the 
fact that in the system there is no frustration (and espe- 
cially no frustration arising from the spatial anisotropics) 

implies that fi = — Z^^'' |/q | — z[^'^ \ga \ , where a ~ x,y ac- 
cording to which axis is characterized by the given mag- 
netic order. 

Without loss of generality, let us consider the situation 
in which the system realizes a parallel antiferromagnetic 
order along the x direction, see Fig. [21(c). All other cases 
can be treated analogously. Following once again the 
steps described in Sec. IIVI and in complete analogy with 
Eq. ([28l) , it is possible to derive the set of conditions that 
the orientation of the E-SQUO must satisfy simultane- 
ously in order for the system to admit a factorized energy 
eigenstate: 

cos2^^ = (53) 

Jx I Jz 

cosH = ^^LZll (54) 

gx - gz 

For the two conditions to be satisfied simultaneously the 
two r.h.s. must coincide and must assume values in the 
interval [0,1]. Hence, by equating them it is possible 
to derive an equation that plays the same role of Eqs. 
(|35]-[3n]). In the hypothesis that the model under investi- 
gation satisfies both Eqs. ([55HM)) . what is left to prove is 
that the associated factorized state be indeed the ground 
state. Following the route illustrated in Sec. IIVI it is 
not difRcult to verify that the following inequalities must 
hold: 

ifx+fy){fx+fz)>0; fx>~fz; (55) 

{gx + gy){gx - gz)>0; gx<gz- (56) 

It is straightforward to establish that they are both sat- 
isfied, provided there are no frustration effects. Namely, 
for spins interacting via couplings of the first kind fa (re- 
ferring to Fig. [21 these arc the intra-chain interactions) , 
the second of Eqs. ([SS]) yields that the denominator in 
the r.h.s. of Eq. (|551) is non-negative and, therefore, the 



corresponding numerator must be non-negative and not 
exceeding the denominator. Hence, we obtain > fy. 
On the other hand, from the first of Eqs. (|55p it fol- 
lows that fx > —fy which, in turn, implies fx > \fv\- 
The latter relation is in agreement with the hypothe- 
sis of absence of frustration. Furthermore, for pairs of 
neighboring spins interacting with a coupling of the sec- 
ond kind ga (referring to Fig. [21 these are the inter-chain 
interactions), the second of Eqs. (|56)) yields that the de- 
nominator in the r.h.s. of Eq. (j54p must be non-positive. 
Hence, also the corresponding numerator must be non- 
positive and therefore g^. < gy Simultaneously, the first 
of Eqs. ([56)) implies gx < —gy and therefore gx < \gy\, 
again in complete agreement with the hypothesis of ab- 
sence of frustration. Collecting all these results we can 
conclude that in the absence of frustration, if a system 
described by the model Hamiltonian Eq. (|5ip satisfies si- 
multaneously all conditions in Eqs. ([55l - [5^ . there exist 
a factorizing field hf at which the ground state is fully 
factorized and characterized by a parallel antiferromag- 
netic order along a;. The factorizing field is determined 
according to the procedure described in Sec. llVi and it 
is not difficult to verify that it reads: 

hf = l^{jr~jr)ijr^jr)- (57) 

The exact form of the factorized ground state reads 

I*/) = (8)l^fc>, (58) 

fc 

l^fc) = (cos(0/2)| t|)+e^'^^sin(0/2)U|)) 

® (cos(0/2)|Ti)+e'^^sin(0/2)|4)) . (59) 

In Eq. ([58|l . the superscript "u" ("Z") denotes the upper 
(lower) chain of the ladder model considered in Fig.[2Uc). 
We observe that (pk = lAI"" does not depend on the choice 
of the chain, implying that the two generic fc-th spins in 
both legs of the ladder arc in the same state equipped 
with an antiferromagnetic order with respect to the near- 
est neighbor sites along the corresponding chain. Analo- 
gous results hold, with the appropriate trivial modifica- 
tions, when considering the other possible magnetic or- 
ders compatible with factorization (See Eq. ([52]) ). 

IX. CONCLUSIONS AND OUTLOOK 

The present work has been motivated both by the 
need for exact solutions to complex many-body quan- 
tum spin models, in particular by the question "when 
is a mean-field solution exact?" , and by the necessity to 
acquire improved knowledge and control on the struc- 
ture of quantum correlations for potential technological 
applications of spin systems. We have introduced a sim- 
ple and rather powerful all-analytic general method to 
determine the conditions for the existence and the prop- 
erties of fully factorized (fully separable) ground states 
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in translationally invariant quantum spin models with 
general two-body exchange interactions and subject to 
external fields. 

The theory developed in the present paper builds on 
and extends the scheme originally introduced in Rcf. 
[131 ■ It provides a readily useful procedure to estab- 
lish, in terms of the Hamiltonian parameters, whether 
the ground state is entangled or completely factorized 
and, in the latter case, what is the exact form of the 
state, the expression of the ground-state energy and of 
the magnetic observables, and the type of ordered phase 
(magnetic order) compatible with factorization. The task 
is achieved by exploiting tools imported from quantum 
information theory such as the formalism of single-qubit 
unitary operations introduced in Refs. (28l . [29| and the 
corresponding entanglement excitation energies, whose 
property of vanishing if and only if a quantum ground 
state is totally factorized plays a central role in the anal- 
ysis of ground-state separability. Using these tools we 
have completed the program initiated in Ref. [13] by de- 
termining a general set of exact relations and inequalities 
that, for frustration-free systems, allow to establish rig- 
orously the occurrence of all types of possible factorized 
ground states. 

We have applied the analytic method to spin-1/2 mod- 
els with general Heisenberg-like exchange interactions of 
arbitrary range (short, finite, long), either isotropic or 
anisotropic, defined on regular lattices of arbitrary di- 
mensions. The method is insensitive to the size of the 
system and applies rigorously in the thermodynamic limit 
as well as for finite lattices. The method allows to estab- 
lish novel exact, fully factorized ground state solutions of 
generally non exactly solvable models, corresponding to 
nontrivial sets of values of the Hamiltonian parameters. 
Key results include, for instance, the rather straightfor- 
ward determination of the existence and form of factor- 
ized ground states in XYZ-type models with different 
types of short-range and long-range interactions on cu- 
bic lattices, a task that would be of formidable complex- 
ity if tackled with numerical techniques or resorting to 
the simple-minded direct method based on the product 
ansatz. 

The complete factorization diagram for nearest- 
neighbor Heisenberg-like anisotropic models might be 
particularly useful for those technological implementa- 
tions which employ spin systems as information proces- 
sors exploiting their ground state entanglement [ij . We 
have provided in Fig. [T] a "minefield" map of the unen- 
tangled working points that need to be avoided, when en- 
gineering the couplings and tuning the external fields, in 
order to achieve satisfactory quantum transmission per- 
formances. On the other hand, there are also alterna- 
tive tasks in quantum information which instead require 
as initial working points factorized states, e.g. for engi- 
neering, via suitable dynamics, specific classes of long- 
distance entangled states useful for quantum state trans- 
fer and quantum tclcportation [4l| or the instantaneous 
creation of strongly entangled graphs or cluster states in- 



volving a mesoscopic or macroscopic number of spins for 
one-way quantum computation [42| : In these instances 
the minefield may turn into a treasure map. The anal- 
ysis of ground-state factorization would be particularly 
useful in the study of quantum spin models on open-end 
lattices. Such systems are not translationally invariant, 
and therefore the analytic method requires to be suitably 
extended to this type of instances. This generalization is 
under way and should be in reach in the near future. 

Notably, according to a general theorem by Kurmann, 
Thomas, and Miiller on factorization [13|, given any 
Hamiltonian of the form Eq. ^ involving higher spins 
S > 1/2, the ground state of that spin-S' system is fully 
factorized at the same value of the external field h ^ hf 
Eq. (j30p , at which factorization occurs in the correspond- 
ing spin-1/2 model. Therefore, the method and the re- 
sults derived in the present paper have a much wider 
scope of application and can be straightforwardly gener- 
alized to interacting systems with arbitrarily high values 
of the spin, provided they arc of the same Hamiltonian 
structure as in the spin-1/2 case. 

Absence of frustration effects, either of geometric or 
dynamical nature, has been essential to the theoreti- 
cal scheme derived in the present work. For instance, 
we have determined the factorization diagram for spe- 
cific examples of frustration-free systems endowed with 
isotropic interaction terms of different spatial ranges, and 
for systems with nearest-neighbor anisotropic interac- 
tions. However, it is clear that there are many extremely 
interesting instances of more complex systems that would 
in general be subject to frustration effects. Continuing a 
stcp-by-step strategy of applying the general method to 
models of increasing complexity, natural further stages of 
investigation would include models subject to frustration, 
e.g., models on linear chains with competing isotropic 
interactions of different range such as a 1-D XYZ-type 
model with nearest-neighbor ferromagnetic interactions 
and next-to-nearcst neighbor antiferromagnetic interac- 
tions. Increasing further the degree of complexity, one 
could then consider models with competing anisotropic 
interactions of different ranges, that according to the type 
of anisotropy could be subject to both geometrically- and 
dynamically- induced frustrations. 

When considering frustrated models, except for very 
special cases, even when the considered system admits a 
fully factorized state as an energy eigenstate, the set of 
conditions derived in the present work are not sufficient 
to establish that such an eigenstate is indeed the ground 
state of the system. In fact, a competition arises be- 
tween factorization, which requires a magnetic ordering, 
and frustration, which tends to destroy magnetic orders, 
ultimately leading to a chaotic correlating behavior. The 
latter unavoidably tends to enhance ground-state entan- 
glement, thus suppressing the occurrence of situations 
where the most favorable state, in terms of energy con- 
tent, is completely uncorrelated. In such cases, a naive 
application of the direct method based on the factorized 
ansatz, for instance to models including anisotropic in- 
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teraction terms of different magnetic types on many dif- 
ferent spatial scales (or even on all scales), leads unavoid- 
ably to overestimate the factorizing effect of the balanc- 
ing between interactions and external fields and thus to 
badly incorrect predictions on the occurrence of fully fac- 
torized ground states in these classes of models, as done, 
unfortunately, in the second part of Ref . [s^] • In a forth- 
coming work (36j we will present a systematic and thor- 
ough study of some important classes of frustrated quan- 
tum spin models in order to establish the conditions for 
the occurrence of true ground state factorization and to 
characterize it in terms of compatibility thresholds with 
frustration. 

An intriguing issue that arises in connection with frus- 
tration is the relationship between factorization and the 
existence of ordered phases. For frustration-free Hamilto- 
nian systems that conserve the even-odd parity, sponta- 
neous symmetry breaking and the existence of an ordered 
phase are necessary conditions for ground state factor- 
ization, because factorized states have no definite parity. 
However, in the presence of frustration there may exist 
factorized energy eigenstates without ground state fac- 
torization. It would then be very interesting to under- 
stand whether in these situations there is still a causal 
relation between factorization and the existence of phase 
transitions in the system. Factorization might then be 
used as an heuristic tool to gain insight in the phase di- 
agram of frustrated quantum systems. 

It would also be important to understand whether chi- 
ral interactions can enhance factorization beyond the lim- 
its imposed by the presence of frustration effects. In fact, 
since factorization can be seen as "mean field becoming 
exact", it seems to require quite naturally a balancing 
between interactions and external fields as the only pos- 
sibility for its occurrence. However, one cannot exclude, 
in principle, that other mechanisms might lead to the 
same effect without requiring the presence of external 
driving fields. 

Future investigations will be concerned with the ap- 
plication of the analytic method to the determination of 
the spectrum of factorization, that is the study of the oc- 
currence of factorized excited states in different classes of 
quantum spin Hamiltonians, and the implications of fac- 
torization diagrams on the structural and informational 
properties of quantum spin models for potential applica- 
tions in quantum technology. 
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APPENDIX: DETERMINATION OF THE 
MAGNETIC ORDERINGS COMPATIBLE WITH 
GROUND-STATE FACTORIZATION 

In Sec. IIVI we have determined the general conditions 
for ground state factorization by selecting as candidate 



factorized ground states only those characterized by a 
ferromagnetic or antiferromagnetic order along the x or y 
directions depending on the value of the minimum value 
/i in the set of the net interactions [see Eq. P^ ]. and 
neglecting all other possible magnetic arrangements. In 
the following, we will prove that indeed, given any spin 
Hamiltonian of the type Eq. ([H]), if a factorized ground 
state exists, then it is characterized by one of the four 
above-mentioned magnetic orders. 

Let us consider the most general factorized state com- 
patible with the constraint of a site-independent magne- 
tization Mz'- 

I*/) = (8) IV'fc) . I^^> = cos(0/2)| U)+e'^i^sin(0/2)| j^) , 
k 

where the local spin state jV'fc) is eigcnstate of 

Ok = cos eSl + cos Lpk sin 6 SI + sin ipk sin OSl (60) 

with eigenvalue equal to 1/2. Let us suppose that there 
exists an Hamiltonian of the type Eq. ^ that admits, 
for some particular value of the external field h = hf, the 
state in Eq. (|5D|) as an eigenstate. Under these hypothe- 
ses we must have that the projection of the factorized 
eigenstate onto the four-dimensional Hilbert space of a 
pair of spins in the lattice must be an eigenstate of the 
pair-Hamiltonian term Hi ^ defined as in Eq. (j23p . 

Following the same strategy as in Sec. IIVI let us intro- 
duce for each site i a set of three auxiliary mutually or- 
thogonal spin operators among which one, that we name 

, coincides with Oi : 

Af ~ cos 9 cos (fiiSf + cos 9 sin (pt^Sf — sin 05f ; 

Af = -sin(^^S'f -f cos(^^5f; (61) 

Al = sin 9 cos '^St + sin 9 sin Lp^S\ + cos 9Sl . 

By inverting Eqs. (j6ip . we can conveniently re-express 
the conventional single-spin operators as functions of the 
new set of operators {Af }: 

Sf = At cos 9 cos (/3j^ — Aj' sin Lp\_ -\- A\ sin 9 cos ^pi, 
Sf = A|cos6'sin(^i + Af cos v3;^-fA| sin 6* sin (/3i^;(62) 
SI = -Af sin 61 + A| cos 61. 

Substituting Eq. (j62p into the pair Hamiltonian we obtain 
the expression for Hi j at factorization as a function of 
the sets of operators {Af} and {A"}, which reads 
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- 2/i^ sin( 



AUA^-- 



cos + sin 9 l^J^. cos cos tpj + Jy sin ipi sin (^j 

+AfA^ [Jl - cos 61) + Af sin sin + cos cos (^^^ - /i^' ^os 6* + Aj^ 

+ 1^ cos 6 + ^1 sin 6^ (^—Jx cos sin tpj + sin cos 
+ (^A^Aj cos 9 + A^A^ sin 0^ ^— sin yjj cos (/Sj + Jy cos sin cpj 



(63) 



where, as usual, r is the distance between sites i and j, 
are the spin-spin couphngs at distance r and along 
direction a, and 

ft-f obeys to the following relation 



2/i} = cos 6I(J: 



cos ifi COS (^j — sin ipi sin 93 j j 



(64) 

^From Eq. (|63p . one has that to ensure that the pro- 
jection of the factorized state onto the four-dimensional 
Hilbert space of spins Si and Sj is an eigenstate of the 
pair Hamiltonian, it is necessary to require that the 
terms associated to the operators AfA^, A^A^, AfA^- 

and A^Af vanish. This condition yields that 
_ — ^ 

cos (fii sin ipj = Jy sin ipi cos ipj , (65) 
J^ sin ipi cos ifj = Jy cos ipt sin ipj . (66) 



Solving Eq. §5^) and Eq. we obtain that if |J^| 7^ 
\Jy\, then the phases f obey to the following relations: 



0,7r/2 ipj ^ ipi + : 



(67) 



with n integer. It is immediate to verify that these rela- 
tions are consistent only with a ferromagnetic or with an 
antiferromagnetic order along the x or the y direction. 
This concludes the proof. 
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